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Abstract — A novel decoding algorithm is developed for general 
quantum convolutional codes. Exploiting useful ideas from clas- 
sical coding theory, the new decoder introduces two innovations 
that drastically reduce the decoding complexity compared to the 
existing quantum Viterbi decoder. First, the new decoder uses 
an efficient linear-circuits-based mechanism to map a syndrome 
to a candidate vector, whereas the existing algorithm relies on a 
non-trivial lookup table. Second, the new algorithm is cleverly 
engineered such that only one run of the Viterbi algorithm 
suffices to locate the most-likely error pattern, whereas the 
existing algorithm must run the Viterbi algorithm many times. 
The efficiency of the proposed algorithm allows us to simulate 
and present the first performance curve of a general quantum 
convolutional code. 

I. Introduction 

We consider quantum error correction coding (QECC) mod- 
els where quantum bits (qubits) and quantum errors are both 
modeled as linear combination of Pauli matrices I, X, Y, Z. An 
[n, k] (block) stabilizer cod^U is defined as a codeword space 6 
which is non-trivially stabilized by a subgroup § of G„, where 
G„ is length-n tensor product of Pauli-operators I, X, Y, Z. By 
definition, the independent generators of § commute to each 
other, and -I i S. 

The dichotomy between the unobservable state of a qubit 
and the observations we can make, lies at the heart of quantum 
computation including QECC. In classic error correction, out- 
put from the channel is observed, measured and used directly 
in the decoding process. Observation in quantum mechanics 
not only destroys the quantum state and makes restoration 
impossible, but the measurement does not tell the original 
quantum state either. The measured outcome is always one 
of the two basic states, whose probability of occurring is 
proportional to the power projection of the original state on 
these two bases. For this reason, to decode a stabilizer code 
typically takes three steps: measuring the syndrome, identify- 
ing the error pattern based on the syndrome, and applying 
a corrective operation to reverse the error. The syndrome 
diagnosis and the error reversal are well-established quantum 
mechanical procedures, and operate the same way regardless 
of the stabilizer configuration. Identifying error patterns, a 
process analogous to syndrome-decoding of classical codes, 
is most challenging. Not only are the techniques used here 
specific to the code structure, but they generally involve 
tedious (exhaustive) search in the Hilbert space. 

Compared to stabilizer block codes, stabilizer convolutional 
codes enjoy the flexibility in codeword length and online 
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We use [...] to denote quantum stabilizer codes, (...) to denote classical 
binary codes, and (...) q to denote classical quaternary codes. 



encoding. The first comprehensive description of the general 
formalism of quantum convolutional codes appeared in (9), 
and a variety of quantum convolutional codes has since 
been proposed JfO- ifTTl . These studies focus on the code 
construction (e.g. based on classical Reed-Solomon codes, 
classical low-density parity check codes, and concatenation 
of classical convolutional codes), and do not usually discuss 
feasible decoding algorithms. The first decoder for quan- 
tum convolutional codes appeared in (9), which developed 
a quantum Viterbi algorithm (QVA) close in spirit to the 
classical Viterbi algorithm. Simplification is made to this 
QVA by attacking only a single unitary error per block (and 
ignoring all the other possible error events) iflOl . IfTTl . For 
quantum convolutional codes consttucted from quasi-cyclic 
sparse matrices, a different type of decoder based on pipeline 
message passing is developed Since the decoder in IflOl . 
full tackles only single errors, and that in |4( performs well 
only on sparse convolutional codes, the QVA in |9) is by 
far the only general and ML decoding method available for 
stabilizer convolutional codes. However, it requires a very 
high complexity, caused by a very large lookup table and 
many rounds of the conventional Viterbi algorithm. The high 
complexity and the induced long delay make this algorithm 
very challenging to implement (and it is for this practicality 
issue that IflOl . fPU proposed to simplify it at the cost of 
degraded performance). 

The contribution of this paper is the development of a 
new decoding algorithm which works for a general quantum 
convolutional code and which is drastically simpler than the 
existing one. Not using lookup tables of any kind, the proposed 
algorithm exploits a simple linear-sequential-circuits based 
mechanism to map a length-(« - k) (binary) syndrome to a 
length-77 (quaternary) candidate vector. This candidate vector 
is subsequently fed to a conventional Viterbi decoder to help 
identify the mostly likely error pattern. The entire process is 
so cleverly engineered that a single candidate vector - any 
one among the pool of 2 (,!+t) possible candidate vectors - 
suffices to locate the error pattern in the ML sense. The high 
efficiency of this algorithm enables us to, for the first time 
in quantum coding literature, simulate and present qubit error 
rate performance curves for a general quantum convolutional 
code. 

II. Quantum Convolutional Codes 

A. Definition and Representations 

A stabilizer convolutional code may be viewed as an infinite 
version of a stabilizer block codes with repeated snucture. The 
stabilizer group, S, for an [n, k, m] stabilizer convolutional code 



is given 



by 0: 

S = sp{M u 



i M ,j, 1 < / < n - k, < ;'}, 



(1) 



where m is the so-called memory parameter, Mqj e G n+mn = 
^/7{I,X,Y,Z} ,8( " + "" ,, . 0, and Mj/s are required to be inde- 
pendent and to commute with each other. Like classical 
convolutional codes, quantum convolutional codes also have 
flexible and adjustable codeword lengths. 

The structure of the stabilizer group generators can also be 
characterized by a semi-infinite matrix M. Similar to the stabi- 
lizer block codes, each line in M represents a generator (some 
Mp), and each column represents a qubit. As illustrated in 
Fig. III-AI M has a block-band structure, and two neighboring 
blocks of generators overlap by mn qubits, representing the 
"actual memory" of the quantum convolutional code. 
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Fig. 1 

Semi-infinite Matrix M of a quantum convolutional code. 



Similar to stabilizer block codes, we can label the polyno- 
mial generators consisting of Pauli elements by F2 polynomial 
vectors in D-domain. The commutation property of stabilizers 
translates to: 

P(D)Q(-f + Q(-)P(D) T = 0, 



D 



D 



(2) 



where D stands for the delay element, similar to that in 
classical convolutional codes, and (P(D)\Q(D)) denotes the F2 
polynomial vectors of a quantum stabilizer code. 

B. Existing Decoding for Quantum Convolutional Codes 

The quantum Viterbi algorithm discussed in (9| is the only 
known general decoder for quantum convolutional codes. A 
syndrome version of the classical Viterbi decoder, the QVA op- 
erates over Pauli-operators (non-binary input). It starts with a 
given list of syndromes, and searches for the candidate vectors, 
which are analogous to the "received sequence (codeword)" in 
classical Viterbi decoding. The task is accomplished through 
table lookup, where the decoder examines the syndromes block 
by block. Because of the convolutional nature, the output 
of two neighboring syndrome blocks overlap by nm qubits, 
and the table lookup procedure must compare and ensure the 
overlapping part of the candidate vectors match. The result is 
a very long list of candidate vectors all of which correspond 
to the given syndromes. Every one of these candidate vectors 
will then be fed into a classical Viterbi algorithm, and among 

in (9) was written as G n + m , which is either mistaken or a typo. 



them, the most-likely sequence, namely, the error pattern, will 
be identified. The overall complexity is rather prohibitive J9). 

To avoid such a high complexity and still be able to 
peek into the code's performance, ifTTl proposed a simplified 
decoder that tackles only the most frequent (but not necessarily 
the most detrimental) error patterns, namely, no more than a 
single error per block. As such, the algorithm is able to inspect 
only two consecutive blocks at a time, and hence drastically 
reduces the size of the lookup table to only 9 terms, at the 
cost of a compromised decoder performance. 

III. Proposed Decoding Algorithm 

A. Syndrome Decoding in Classical Coding 

The proposed decoding algorithm roots to the classical 
syndrome decoding approach for linear codes. In essence, 
syndrome decoding is minimum distance decoding, which is 
equivalent to the maximum likelihood decoding if the channel 
is discrete memoryless with error probability strictly less than 
0.5. The linearity of a code allows the code space to be 
arranged in a special way, termed the standard array in coding 
jargon. For an («, k, d) classical linear block code, the standard 
array groups all the 2" n-bit vectors into 2 k columns and 2"~ k 
rows, such that (i) the vectors in the same row corresponding to 
the same syndrome, (ii) the first row consists of the set of valid 
codewords (starting with the all-zero codeword), and (iii) the 
first column consists of all the coset leaders, i.e. correctable 
error vectors with Hamming weights lower than \_{d - 1)/2J. 
Each row of the standard array is called a coset. In concept, 
the function of a syndrome decoder is to identify the coset 
leader, i.e. the minimum-weight vector, corresponding to a 
given syndrome, but this seemingly simple function can be 
extremely expensive to implement in practice. 

B. Decoding Algorithm 

Stabilizer codes in quantum coding are analogous to lin- 
ear codes in classical coding, and may be tackled through 
a similar coset representation and syndrome decoding. The 
propose decoder here is not only a conceptual quantum version 
of the syndrome decoder, but also a practical mechanism 
that provides a systematic and concrete way and practical 
procedure to realize the syndrome decoder with manageable 
complexity. The proposed quantum syndrome decoding, which 
is a low-complexity ML decoding procedure, is depicted in 
Fig. |2] Below we explain why, how, and how efficiently 
the proposed algorithm works. We start by transforming the 
stabilizer generators to a binary polynomial form, and find its 
equivalent transfer polynomial. We then detail the step-by-step 
approach in Fig. [2] and demonstrate examples to show how to 
decode without (exhaustive) lookup tables. 

Here are a few notations used in the discussion. A quantum 
convolutional code may be described by its semi-infinite sta- 
bilizer matrix M, or stabilizer polynomial S(D). The proposed 
decoding procedure involves the derivation of the equivalent 
F2 or F4 classical code from the stabilizer polynomial S(D), 
which are termed the equivalent ¥2 or F4) transfer polynomial, 
and denoted by Hb(D) or H q {D), respectively. The decoding 
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Diagram of general quantum decoding, composed of two stages: decoding error patterns, recovery. 



procedure also makes heavy use of useful concepts in classical 
coding: syndrome former (SF) and inverse syndrome former 
(ISF). For a given (n,k) linear code and hence a deterministic 
standard array, the SF finds the syndrome associated with a 
given sequence (an error pattern), and the ISF finds one error 
pattern - any one of the 2 k possibilities - that associates with 
a given syndrome. 

Decoding Algorithm: 

The "syndrome decoder" block (dashed box) in Fig. ^illus- 
trates the proposed quantum syndrome decoder, whose input 
is a binary syndrome sequence obtained from the syndrome 
measurement. 

1) Represent the original quantum convolutional code using 
its equivalent F2 transfer polynomial, which leads to an 
equivalent classical convolutional code. 

2) Perform syndrome decoding on the equivalent classical 
convolutional code, through the following two substeps: 

a) Compute the F2 ISF from the equivalent F2 transfer 
polynomial. Use the ISF to locate an arbitrary 
candidate vector in the coset corresponding to the 
given syndrome. 

b) Compute the F2 generator polynomial Gb(D) from 
the equivalent F2 transfer polynomial. From the 
generator polynomial, implement a conventional 
trellis decoder (such as the Viterbi algorithm or 
the BCJR algorithm), and use it to decode the 
candidate vector into a valid codeword. Exclusive- 
OR (XOR) the candidate vector and the valid 
codeword to get the error pattern. 

The output of the decoder is a 277-bit binary sequence, which 
can be re-assembled to an n-qubit Pauli-formed error pattern. 
This error pattern is then fed into the "recovery" block in Fig. 
H] and counter-acted to recover the original quantum states. 

Step 1: Deriving equivalent F 2 transfer polynomial 

Let S (D) be the stabilizer polynomial of a quantum convo- 
lutional code. Recall that a stabilizer matrix can be described 
either in the binary F2 form, S = (P\Q), where P and Q are 
symplectically orthogonal, or in a quaternary F4 form, where 
Pauli operators {I, Y, X, Z} are represented by F4 elements 
{0, 1, co, co 2 = to}, respectively, conforming to the zero trace 
inner product (to is a primary element in F4). A stabilizer code 
in F4 form takes a similar flavor as an F4 classic code, but the 
syndrome takes a binary, rather than quaternary, value. Let 



Mj - (m\j,m2j) be a binary row in (P\Q), and let E — (e\,e%) 
be a binary error vector. Without loss of generality, re-write 
E and My in F4 as: E — e\ + eito, Mj - my + m%jto. The 
syndrome corresponding to E is computed as 

Sj(E) = e\ ■ m\ + g2 ■ m\ ■ (3) 

An immediate and important implication of (01 is 

Lemma 1: J4) An equivalent classical linear code to the 
quantum convolutional code with stabilizer generator S(D) = 
(P(D)\Q(D)) has a binary (F 2 ) transfer polynomial: H b (D) = 
P(D 2 ) + DQ(D 2 ). Here equivalent is meant in terms of de- 
coding capability promised by the code (disregarding possible 
degeneracy). 

Step 2: Syndrome decoding for convolutional codes 

The problem now boils down to performing syndrome 
decoding for the equivalent classical convolutional codes. 

For an (n, k) linear code, the function of the SF is to map a 
length-n vector to a length-(n-Af) syndrome. The ISF performs 
the opposite function of SF by locating a length-n vector 
in the coset associated with a given syndrome. Since each 
coset contains 2 k vectors, there are 2* possible outputs an ISF 
may produce, and the resultant vector can be any of the 2 k 
vectors depending on the specific inverse syndrome former 
being implemented. The choice of the ISF will not affect the 
end-result of the proposed syndrome decoder. Efficient and 
systematic ways exist to implement SF and ISF, the simplest 
of which is through linear transformation or linear circuits. 
For a linear block code, the parity check matrix can serve the 
role of SF, and the left inverse of the parity check matrix can 
serve as a matching ISF. 

The proposed syndrome decoder makes essential use of 
the functions of SF and ISF. The "syndrome measurement" 
block in Fig. |2] generates a measure of the syndrome, and the 
"syndrome decoder" block uses a matching ISF to acquire a 
valid vector in the coset associated with that syndrome, and 
subsequently locates the coset leader after trellis decoding 
and an XOR operation. Finally, the odd-positioned bits and 
the even-positioned bit in the coset leader are paired up to 
obtain the tensor product form I, X, Y, Z of the error pattern, 
which will then be applied to the received quantum sequence 
(the "recovery" block in Figj2] a well-established quantum 
mechanics procedure) to get the correct codeword. 

Lemma 2: The proposed syndrome decoding approach in 
Figf2]is valid. 



Proof: The validity of this decoder is warranted by the fact that 
the error pattern/coset leader is the minimum-weight vector 
inside each coset. Suppose that the transmitter sends a valid 
binary codeword U. The discrete memoryless channel adds 
some recoverable binary noise e to U to yield V. That is, 
V = U ffi e, and hence e — U ffi V, where © denotes binary 
addition (XOR). The quantum receiver collects V, but instead 
of having a direct measurement of V, it gets a measurement 
of the syndrome S . We show the the combination of ISF and 
trellis decoding in Fig. [2] successfully deduce e. 

Inverse syndrome former maps S to an arbitrary vector W 
in that coset. In the standard array, U locates in the first row 
(corresponding to the zero syndrome, V, e and W all locate in 
the same row corresponding to the syndrome S , e is the coset 
leader (because it is a recoverable error patten), and V locates 
in the same column as U (because E/ffie = V). Hence, W is in 
fact some codeword corrupted by error e, namely, there exist 
some valid codeword N such that N ffi e = W. 

The vector W at the output of the ISF is fed to the ML 
trellis decoder. If the code can support this channel, then with a 
probability approaching 1, the trellis decoder is able to deduce 
the valid codeword N from W. Since e = U ffi V = N ffi W, 
XORing N (output from trellis decoder) and W (output from 
ISF) thus produces the target error pattern e. □ 

C. Realization of Important Modules in Decoder 

As discussed before, the proposed syndrome decoder in 
Fig. [2] involves two key modules, the ISF and a conventional 
trellis decoder for the equivalent classical convolutional code. 
The latter is a well-established procedure, given a generator 
polynomial. Lemma 1 establishes the relation between a 
stabilizer polynomial S(D) of a quantum convolutional code 
and the transfer polynomial H/,(D) of its equivalent classical 
linear convolutional code, below we focus the discussion on 
the derivation of the ISF and the generator polynomial for a 
given H h (D). 

There exists intrinsic connection between F2 generator poly- 
nomial, transfer polynomial, SF and ISF. Specifically, consider 
a binary k/n convolutional code with transfer polynomial 
Hi,(D). The transpose of the transfer polynomial, (Hb) T , which 
has dimension « x (« - k) can serve as an SF. The generator 
polynomial can be implemented using a linear sequential 
circuit specified by k x n transfer function Gb with rank k 
satisfying 

G b x (H b ) T = k , (4) 

where 0^ is a k-by-(n - k) all-zero matrix. The constraint in (0]i 
guarantees that all the valid codewords are associated with the 
length-(« - k) all-zero syndrome 0,,-t and that a set of length-n 
codewords/vectors have the same syndrome if and only if they 
belong to the same coset. 

Inverse syndrome former can be obtained by taking the left 
inverse of the syndrome former, and is therefore denoted by 

i»b l ) T -- 



l n -k, 



(5) 



where is an identity matrix with rank n-k. The inverse 
syndrome former may be derived using the partial Gaussian 



elimination or the Moore-Penrose matrix inverse method. The 
latter states that the left inverse of a rectangular matrix B, 
denoted as B + , takes the form of B + = (A 7 -A) -1 -A T . Although 
the Moore-Penrose matrix inverse was originally derived in 
real-valued matrix, it can be easily extended to binary matrix 
by with module-2 computation. 

Example 1: This example illustrates how to obtain the ISF 
and the generator polynomials from a given SF. Supposed the 
given SF takes of the form of 



l+D 2 
l+D + D 2 



ISF 
GP 



which corresponds to a rate 1/2 (classical) recursive convolu- 
tional code. Using the Moore-Penrose procedure, a possible 
ISF and generator polynomial may be obtained: 

(h;V = ( i + d D); 

For a given SF, matching ISFs and generator polynomials 
are not unique. For instance, in the previous example, another 
possible choice for ISF and generator polynomial are: 



(6) 
(7) 



ISF: 
GP : 



(HbV 



I 1 1+Q 2 2 ); 
V l+D+D 2 J ' 



Gb = ( l+D + D 2 l+D 2 Y 



(8) 
(9) 



It is desirable to choose a non-catastrophic generator polyno- 
mial, and an inverse syndrome former with a low complexity. 

Example 2: This example illustrates the entire decoding 
procedure. Consider the quantum convolutional code in ifTTIl . 
which has a semi-infinite stabilizer matrix (in Pauli operator 
form): 
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(10) 



It represents a [n,k,m] = [3, 1, 1] quantum convolutional code 
with rate 1/3. 

Following the proposed syndrome decoding algorithm, the 
first step is to derive the equivalent F2 transfer polynomials: 



H h (D) 



l+D 2 



1 +D' 
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D + D 2 



(11) 



whose matrix transpose, Hb{D) T , represents the syndrome 
former of a (3, 1,2) F2 (classical) convolutional code. 

The next step is to derive the ISF (H^(D)y l and the 
generator polynomial Gi,(D) matched to the given SF Hb(D) T 
(all in F2). Following the Moore-Penrose procedure, we get 

_ 1 L 

ISF : 



(HliD))- 1 



D+D 
1 







D 2 +D' D : +D'+D 



D+D 2 +D 3 





(12) 



GP: G b (D) = [D 2 l+D 2 l+D 2 ). (13) 

The resultant decoder diagram is shown in Fig. [3] The 
binary syndrome sequence S 1 is separated into two streams, 



and feeds into the inverse syndrome former, which, in this 
example, is like a recursive convolutional encoder with two 
streams in and three streams out. The output streams assemble 
to a sequence W, passes through a BCJR decoder, and gen- 
erates a sequence which is a valid codeword. The difference 
between the codeword sequence and W is the error pattern 
sequence e. e is a binary sequence, and will be divided into 2 
parts and packed back to Pauli operators accordingly. 
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Fig. 3 

Practical decoder of quantum convolutional code (example 

2). 



D. Decoder Implementation in ¥4 

The proposed decoder is applicable to general quantum 
convolutional codes on a variety of quantum channel models 
characterized by memoryless independent Pauli operators. The 
previous discussion and examples are centered around F2 
implementation. When the quantum convolutional code is 
derived from an F/t-linear codes, such as those proposed in 
ifTTI . one may either find its equivalent binary convolutional 
code as discussed before, or explore syndrome decoding on its 
equivalent F4 counterpart. The latter is possible, because the 
syndromes, which are typically («-A:)-bit binary sequences, 
may also be represented by a length (n-k)/2 ¥4 sequences in 
this example. 

Example 3: Consider the rate- 1/3 [3, 1, 1] quantum convolu- 
tional code from Example 2, whose equivalent F4 convolu- 
tional code has a F4 transfer polynomial H q {D) 



HJD) = (1 + D, 1+cjD, 1 + uD). 



(14) 



The F4 syndrome former is given by H q {D) T . 

The ISF and the generator polynomial may also be derived 
in F4 form: 



ISF: (H q \D)) T = ( 1, 1, 1 ) 
GP : 



(15) 



. 0, 1 +a>D, 1 +u)D 
G * (jD) = ll + w D, 1+A l- (16) 



Clearly, the trellis decoder (Viterbi/BCJR) in Fig. |2 can be 
implemented based on F4 generator polynomial. The entire 
syndrome decoding is carried out over a (n,k,m) = (3,2,1) 
classical F4 convolutional code. 



IV. Simulation Results 

We demonstrate the proposed algorithm by simulation on a 
memoryless bi-polar Pauli error channel, where X and Z occur 
at probability p - p 2 , and Y occurs at probability p 2 . 



Pr(X) = Pr(Z) = p-p , Pr(Y) = p . add 6 "I" padding after each frame 




Fig. 4 

Simulation results of quantum convolutional codes. 

Notice that the ISF is essentially a recursive convolutional 
code, padding bits are necessary to generate correct W for 
simulations. For an [n, k, m] quantum convolutional code, the 
padding sequence should be a length-(« + nm) I sequence. We 
simulated the [3,1,1] quantum convolutional code at a frame 
length of 900 qubits and tested 10000 frames. This is like 
a [900, 300] stabilizer block with an additional 6-bit padding 
(overhead), so the exact code rate is 300/(900 + 6) « 1/3. The 
resultant performance curve is shown in Fig. @] 

V. Conclusion 

We have developed a new syndrome decoding approach, 
which is applicable to a general stabilizer convolutional code 
on a general Pauli-error channel model. Comparing with the 
existing quantum Viterbi decoder and its simplified version, 
the new algorithm enjoys significantly lower complexity (no 
lookup table, and one-time Viterbi) without scarifying ML 
performance. The decoding algorithm roots to the classical 
syndrome decoding theory, and boasts provenly optimal per- 
formance. Using this algorithm, we simulated and presented 
the first performance curve for a general quantum convolu- 
tional code. 
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